A NOTE ON SHELLING 
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Abstract. The radial distribution function is a characteristic geometric quantity of a point 
set in Euclidean space that reflects itself in the corresponding diffraction spectrum and 
related objects of physical interest. The underlying combinatorial and algebraic structure is 
well understood for crystals, but less so for non-periodic arrangements such as mathematical 
quasicrystals or model sets. In this note, we summarise several aspects of central versus 
averaged shelling, illustrate the difference with explicit examples, and discuss the obstacles 
that emerge with aperiodic order. 



1. Introduction 

One characteristic geometric feature of a discrete point set A C M d , which might be thought 
of as the set of atomic positions of a solid, say, is the number of points of A on shells of radius 
r around an arbitrary, but fixed centre in Of particular interest are special centres, such 
as points of A itself, or other points that are fixed under non-trivial symmetries of A. This 
leads to the so-called shelling structure of A. Here, we consider infinite point sets only. In 
general, one obtains different answers for different centres, and one is then also interested in 
the average over all points of A as centres, called the averaged shelling. 

The spherical shelling of lattices and crystallographic point sets (i.e., periodic point sets 
whose periods span ambient space) is well studied, and many results are known in terms of 
generating functions. If A is a lattice, the number of points on spheres of radius r centre 
(central shelling) is usually encapsulated in terms of the lattice theta function |16[ Ch. 2.2.3] 

(i) e A (z) = £y = 

x£A k 

where q = e mz and c(k) is the number of lattice points of Euclidean square norm (= square 
length) k. A closed expression for the latter can be given in many cases, see |161 Ch. 4] 
for details on root and weight lattices and various related packings, and [Hj for an explicit 
example. There are many related lattice point problems, see |24j and references therein for 
recent developments. 

One special feature of a lattice is that the shelling generating function is independent of the 
lattice point which is chosen as the centre - and, consequently, central and averaged shelling 
give the same result. Similarly, for uniformly discrete point sets that are crystallographic, 
the average is only over finitely many points in a fundamental domain and can often be 
calculated explicitly. For general uniformly discrete point sets, however, the situation is more 
complicated in that no two centres might give the same shelling function, or that the average 
may not be well defined. But there is one important class of point sets, the so-called model sets 
(also called cut-and-project sets, see (HD EJ EH HH1 EEl EE31 EES] and references therein), which 
provide a high degree of order and coherence so that an extension of the shelling problem 
to these cases is possible, and has indeed been pursued. The original motivation for the 
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investigation of model sets came from applications in physics. Meanwhile, due to interesting 
connections with several branches of mathematics, they are also studied in their own right, 
see |3f)[ I37( 111 j and for details and further references. Below, we shall summarise the key 
properties of model sets needed for this article. 

One of the earliest attempts to the shelling of model sets, to our knowledge, is that of Sadoc 
and Mosseri [30] w h° investigated the 4D Elser-Sloane quasicrystal and then conjectured 
a formula for the central shelling of a close relative of it which was obtained by replacing the 
highly symmetric 4D polytype used in ^H] by a 4D ball as a window. The conjecture was put 
right and proved in 33,. by means of algebraic number theory revolving around the arithmetic 
of the icosian ring I, a maximal order in the quaternion algebra HI(Q(v / 5 ))• Recently, the 
central shelling was extended to the much more involved 3D case of icosahedral symmetry |46l . 
Also, some results exist on planar cases, e.g., for special eightfold and twelvefold symmetric 
cases with circular windows J33J . 

The common aspect of all these extensions to model sets (or mathematical quasicrystals) is 
that only the central shelling of a highly symmetric representative has been considered, with 
a ball as window in internal space. This is a rather special situation which appears slightly 
artificial in view of the fact that the most relevant and best studied model sets usually have 
polytopes rather than balls as window, or, more generally, even compact sets with fractal 
boundary, cf. 4 . Also mathematically, the classical examples such as the rhombic Penrose or 
the Ammann-Beenker tiling are very attractive due to their rather intricate and unexpected 
topological nature |3| IT9]. 

A more natural approach to model sets seems to be the averaged shelling, and it is the 
aim of this article to start to develop this idea. As we shall see, the topological structure will 
be manifest in the examples discussed below. On the other hand, the central shelling does 
have a universal meaning, too, if one considers it first for modules rather than for model sets. 
The window condition can then be imposed afterwards, see [H1IS1IZI f° r some examples. This 
approach is implicit in [SSj, but does not seem to have attracted much notice. It is important, 
though, because it leads to a separation of universal and non-universal aspects. 



2. Central shelling 

Many of the well studied planar tilings with non-crystallographic symmetries share the 
property that their vertices (or other typical point set representatives) form a discrete subset 
of rings of cyclotomic integers. This gives a nice and powerful link to results of algebraic 
number theory, which has, in fact, been used to construct model sets |38| . and which also 
appeared before in a different context |29| . Let us thus first explain the situation of central 
shelling for these underlying dense point sets. 

Let £ n be a primitive n-th root of unity (with n > 3), e.g., £ ra = e 27 ™/ n , and Q(£ n ) the 
corresponding cyclotomic field. Then, Q(£ n + £ n ) is its maximal real subfield. From now on, 
we will use the following notation 

(2) K = Q(£„), k = Q(£ n +4), = Z[6J, o = Z[£« + e„], 

where O is the ring of cyclotomic integers, which is the maximal order of K, and o is the ring 
of algebraic integers of k, see |35j . 
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Note that O is a Z-module of rank (f(n), where (p denotes Euler's totient function. The set 
O, seen as a (generally dense) point set in M 2 , has iV-fold rotational symmetry, where 



(3) N = N(n) 



if n is even, 
if n is odd. 



This also means that O has precisely N units on the unit circle, which are actually all roots of 
unity of K. Also, IK is a totally complex field extension of k of degree 2. It is known that, in 
this cyclotomic situation, the unique prime factorisation property of O (i.e., class number one) 
implies that of o, and this happens in precisely 29 cases, compare [451 Thm. 11.1], namely 
for 

n £ {3, 4, 5, 7, 8, 9, 11, 12, 13, 15, 16, 17, 19, 20, 21, 
^ 24, 25, 27, 28, 32, 33, 35, 36, 40, 44, 45, 48, 60, 84} , 

where n ^ 2 mod 4 to avoid double counting. Note that n = 1 (N = 2) is excluded here 
because it corresponds to IK = Q with O = Z, which is only one-dimensional. 

Now, let p be a prime of o. Then, in going from o to O, precisely one of the following 
cases applies, see (HIJ Ch. I] and [13 Ch. 4]: 

(1) p ramifies, i.e., p = PP with P a prime and P/P a root of unity in O. 

(2) p is inert, i.e., p is also prime in O. 

(3) p is a splitting prime of O/o, i.e., p = PP with P/P not a unit in O. 

Up to units, all primes of O appear this way. 

Prime factorisation in o versus O can now be employed to find the combinatorial structure 
of the shells. We encode this into the central shelling function c(r 2 ) which counts the number 
of points on shells (circles) of radius r. By convention, c(0) = 1. 

Theorem 1. Let O = Z[£ n ] be any of the 29 planar 7L-modules that consist of the integers 
of a cyclotomic field with class number one. Then, for r 2 > 0, the function c(r 2 ) vanishes 
unless r 2 E o and all inert prime factors of r 2 occur with even powers only. In this case, 

(5) c(r 2 ) = N- J] (t(p) + l), 

p\r 2 
p splits 

where p runs through a representative set of the primes of o. Here, t(p) is the maximal power 
t such that p l divides r 2 . The pref actor, N = N(n) of Eq. (J3J), reflects the point symmetry 
of the module. Furthermore, r 2 is then a totally positive number in o, i.e., all its algebraic 
conjugates are positive as well. 

Proof. Since c(0) = 1 by convention, consider r 2 > 0. If there exists a number x <E O on the 
shell of radius r around 0, we must have r 2 = xx, hence r 2 G o. In this case, any inert prime 
factor p of r 2 (in o) necessarily divides both x and x (in O). Consequently, the maximal 
power t = t(p) such that p l divides r 2 must be even. 

Conversely, assume r 2 > and t(p) even for all inert primes of o. If a ramified or a 
splitting prime p = PP divides r 2 , we know that equal powers of P and P occur in the prime 
factorisation of r 2 in O. Consequently, we can group the prime factors of r 2 in O into two 
complex conjugate numbers, i.e., we have at least one solution of the equation r 2 = xx with 
x £ O, so c(r 2 ) > 0. 



4 



MICHAEL BAAKE AND UWE GRIMM 



Consider a non-empty shell with r 2 > 0, i.e., r 2 = xx for some / i 6 O. Consider the 
prime factorisation r 2 = e • pj 1 • . . . • pi in e>, with e a unit. If pj is not a splitting prime, the 
distribution of the corresponding primes in O to x and x is unique, up to units of O. 

If, however, pj = PjPj is a splitting prime, we have to distribute {PjPjfi over x and x. 
In view of Pj being the complex conjugate of Pj, but not an algebraic conjugate, we have the 
options of (Pj) s (Pj) t: >~ s as factor of x and (Pj) s (Pj) t:i ~ s as factor of x, for any < s < tj. 
This amounts to t- + 1 different possibilities, which gives the corresponding factor in ©. 

As mentioned above, there are N units of O on the unit circle. This means that, as soon 
as r 2 > 0, points on the shells come in sets of N, which gives the prefactor in Q. Together 
with the previous arguments, this explains the multiplicative structure of c/N. 

Finally, assume r 2 = xx for some ^ x £ O and let a be any Galois automorphism of IK 
over Q. Then we have 

< a(x)a(x) = a(x)a(x) = a(xx) = a(r 2 ) 
so also all algebraic conjugates of r 2 are positive. This shows that r 2 is totally positive. □ 

Remark: It is clear that Theorem^ can be generalised to the situation that K is a totally 
complex field extension of a totally real field k whenever IK has class number one, with sets of 
integers O and o as above, compare |45l Thm. 4.10]. In this case, the prefactor in Eq. Q has 
to be replaced by the number of elements in the unit group of O that lie on the unit circle. 

Let us consider the cyclotomic case in more detail. If V(r 2 ) = {x £ O \ xx = r 2 }, and 
a is any Galois automorphism of K/Q, then c{0) = O and V(r 2 ) is mapped bijectively to 
V(o(r 2 )). This means that c(r 2 ) = c(a(r 2 )). 

Moreover, consider the situation that two totally positive numbers of o, r 2 and R 2 , are 
related by R 2 = er 2 , with e a unit in o. Clearly, e is then also totally positive. If e is of the 
form e = uu, with u a unit in O, the mapping x i— > ux gives a bijection between V(r 2 ) and 
V{R 2 ), hence c(r 2 ) = c(R 2 ). If all totally positive units of o are of this form, which includes 
the case that e is the square of a unit in O, we may conclude that the central shelling function 
c only depends on the principal ideal of o generated by r 2 . 

In general cyclotomic fields, this factorisation property of totally positive units need not 
be satisfied (e.g., it fails for n = 29). However, it is true for all class number one cases. More 
precisely, if n is a power of 2, all totally positive units of O are squares of units of O, which is 
known as Weber's theorem, compare |211 Cor. 1 and Rem. 2]. The same statement holds if n 
is an odd prime below 100, except for n = 29, see |21| Ex. 2]. We checked explicitly, using the 
KANT program package |17l I26j , that this remains true for all n from our list © that are 
prime powers. All remaining cases of (@} are composite integers. Here, not all totally positive 
units of O are squares in O, but they are of the form e = uu, with u a unit in O. This was 
again checked using KANT. The difference to the other cases comes from the additional unit 
u = 1 — £ n , compare (33 Cor. 4.13]. 

We may conclude as follows. 

Fact 1. Let O = Z[£ n ] be any of the 29 planar "L-modules that consist of the integers of 
a cyclotomic field with class number one, and o = Z[£ n + £„]. Then, the central shelling 
function c for O depends on r 2 G o only via the principal ideal r 2 o generated by it. □ 

This allows us to reformulate the result of Theorem ^ by means of ideals and characters 
of the field extension K/k. By a character \ ^ 0, we here mean a totally multiplicative real 
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function of the ideals of O, i.e., x( a ^) = x( a )x(b) for all ideals o and b of O, see (SHI Ch. VII. 6] 
for background material. In particular, x{°) = 1- It suffices to specify the values of \ f° r an 
prime ideals p of c We define 



(6) X(P) 



if p ramifies 
-1 if p is inert 

1 if p splits 

where the property of the prime ideal p refers to the behaviour under the field extension from 
k to K. This leads to the following result. 

Corollary 1. Under the assumptions of Theorem^ the central shelling function c is propor- 
tional to the summatory function of the character \ of Eq. ©, i.e., 

(7) c(r 2 o) = N- J2 X(fl), 

a\(r 2 o) 

with N given by Eq. J3J. 

Proof. Due to unique prime factorisation in o and the multiplicative structure of c/N accord- 
ing to Eq. (JSJ), it is sufficient to verify the claim for prime powers, i.e., for r 2 o = p e . Clearly, 
if p ramifies, the sum in Eq. (JJJ) gives c(p^) = 1 for all £ > 0. The alternating sign of x(p ) 
for inert p implies c(p l ) = for odd I and c(p^) = 1 otherwise. If p splits, the right hand side 
of Eq. (J7J) adds up to I + 1. Invoking Fact ^ an d a comparison with Eq. (JSJ completes the 
proof. □ 

The explicit use of Theorem ^ and Corollary ^ requires the knowledge of the splitting 
structure of the primes. Examples can be found in j^H El EI, see also [3IJ If one is 
interested in the central shelling of a model set rather than that of the underlying (dense) 
module, one has to take the window into account as a second step. A model set A(JT) in 
"physical space" M. d is defined within the following cut-and-project scheme |31l |1] 



(8) 





R d x H 




H 




> 


U 


U 




U dense 


L < 


r 


> 


L* 



where the "internal space" H is a locally compact Abelian group, and r C R d x H is a 
lattice, i.e., a co-compact discrete subgroup. The projection L* = n H (r) is assumed to be 
dense in internal space, and the projection into physical space has to be one-to-one on r. 
Consequently, the mapping *: L — ► L* C H, with * = -k h o (vr| r ) , is well defined. It is 
called the *-map of the cut-and-project formalism, compare |31j . Note that the *-map need 
not be injective, i.e., its kernel can be a nontrivial subgroup of L. 
A model set A(f2) is now defined as 

(9) A{Q) = {x G L | x* e Q} = {ir{y) \ y G T, 7r H (y) e Q} C R d , 

where the window Q C H is a relatively compact set with non-empty interior. Usually, one 
either takes an open set or a compact set that is the closure of its interior. Note that the 
*-map is well defined on A(f2), with (A(f2)) C fi. More generally, also sets of the form 
t + A(Q) with t G R d are called model sets. If t G L, one has t + A(Q) = A{t* + Q) and is 
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back to the case of Eq. Q, which is sufficient for our discussion. For the above example of a 
cyclotomic field K = Q(£n)> we need d = <p(n) to construct model sets with n-fold symmetry, 
compare ^Ql App. A]. 

In order to compute the central shelling for a model set A(Q), one first determines all 
points of the module L = vr(7~') on the shell of a given radius r. Then, the window fl decides, 
according to the filtering process of Eq. @, which of these points actually appear in the 
model set, and the shelling formula is modified accordingly. As long as we are dealing with a 
one-component model set (i.e., as long as all points are in one translation class), the formula of 
Theorem ^ thus gives an upper bound on the shelling number in the model set. As mentioned 
above, the central shelling of a few model sets with spherical windows |4U1 1331 1341 1351 146j has 
been considered in detail. 



A moment's reflection reveals that the averaged shelling is considerably more involved. In 
order to determine the averages, one would need to know all possible local configurations 
up to a given diameter together with their frequencies, provided the latter are well defined. 
In general, this is not the case, as cluster or patch frequencies in general Delone sets need 
not exist. However, regular model sets are particularly nice in this respect because all patch 
frequencies exist uniformly which is equivalent to unique ergodicity of the corresponding 
dynamical system |44| 142) (under the translation action of M. d ). Moreover, due to existence 
of the cut-and-project scheme (jSJ) and Weyl's theorem, compare |32j . it is possible to transfer 
the averaging part of the combinatorial problem to one of analysis. 

Let us also point out that Eq. for a model set does not make much sense as it would 
depend on the representative chosen, rather than being a quantity attached to an entire 
local indistinguishability (LI) class, compare |431 0]. If A is a lattice, A — A = A, and 
we could equally well sum over the difference set in (Q). Using this for model sets would 
give Y^xeA-A Q X X w hich is constant on the LI class. However, this still does not reflect the 
statistical aspects of the (local) shells, because each x S A — A is counted with weight one. 
Let us thus introduce the averaged shelling function a(r 2 ) as the number of points on a shell 
of radius r, averaged over all points of A as possible centres of the shells. 

Now, let A = A(f2) be a regular, generic model set, in the terminology of [HJ|) with window 
]?, i.e., Q is a relatively compact set in H with non-empty interior, boundary of measure 0, 
and dfi n ir H (r) = 0. For simplicity, we also assume that H = R m , though a generalisation 
of what we say below to the case of general locally compact Abelian groups is possible. In 
analogy to Eq. (Q), a generalised theta series could be defined ad hoc as 



where q = e mz and 7Z = {r £ M>o | \y\ = r for some y £ A — A} is the set of possible radii as 
obtained from the set of difference vectors between points of A. The coefficient a(r 2 ) is now 
meant as the averaged quantity defined above, which we will now calculate. 

Let v(y) denote the relative frequency of the difference y between two points of the model 
set (hence y £ A — A). Up to the overall density of the model set, v{y) is an autocorrelation 
coefficient of the point set A. This quantity exists uniformly for all y as a consequence of the 
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model set structure (ISllfflEU- But then, we obviously obtain 

(11) a(r 2 ) = £ u(y). 

yeA-A 
\y\=r 

On the other hand, if A s = {x £ A \ \x\ < s}, one has 

u(y) = lim — — r 1 = lim , . , . , 1 

yy> 8 -0O MJ ^ s->oo \ (A S )*\ ^ 

x+yGA ( x+y y £ Q 



(12) = -rr^ / l n (z)l a (z + y*)dz 

where 1^ is the characteristic function of the window. Note that, as the *-map need not be 
injective, the second equality may only hold in the limit s — > oo (this step is implicit in the 
proof of |41l Thm. 1]). We add it here because it shows how the counting is transfered to 
internal space, in particular in the cases where the *-map is one-to-one, which is the situation 
we will meet in the examples. 

The last step in (|12|) is now a direct application of Weyl's theorem on uniform distribution. 
This is justified here because (A s )* , for increasing s, gives a sequence of points in Q that 
are uniformly distributed, see [23 E2j and [12 Thm. 1], and because dQ has measure by 
assumption. In this situation, the averaged quantities are the same for generic and singular 
members of the LI class |411 Ej. Moreover, it also does not change if dQ n L* ^ 0, so that 
the corresponding assumption can be dropped. Consequently, the averaged shelling function 
is constant on LI classes of regular model sets. We combine Eqs. (jllj) and (|12[) to obtain 

Theorem 2. Let A be a regular model set in the sense of Moody |31| . obtained from a cut- 
and-project scheme © with internal space H = R m and window Q. Then, the averaged 
shelling function a(r 2 ) exists, and is given by 

(13) a(r 2 ) = - vol(/2 n(/2 -</*)). 

^ ' y^A-A 
\y\=r 

In particular, a(r 2 ) vanishes if there is no y G A — A with y ■ y = r 2 . □ 

Remark: This result allows the calculation of the shelling function, for any possible radius 
r, by evaluating finitely many volumes in internal space. This is so because a model set A 
has the additional property that also its difference set, A — A, is uniformly discrete, so that 
there are only finitely many different solutions of \y\ = r with y E A — A. 

4. Examples 

Let us first consider a well-known model set in one dimension, the Fibonacci chain, which 
can be described as 

(14) A F = {x£Z[r] | x* G [_1, T -1]} = A([-1,t-1]), 

where Z[t] = {m + nr | m,n £ Z} is the ring of integers in the quadratic field Q(t) and 
t = (l + \/5 )/2 is the golden ratio. The *-map in this setting is algebraic conjugation in Q(r), 
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defined by y/E i— > —y/E. The 2D lattice behind this formulation is r = {(x, x*) \ x € Z[r]}. A 
short calculation results in Ap — A-p = A{\\— t, t]), and 



(15) v(y) = v{-y) = f F (y*) 



if \y*\ > t 
It if \y*\ < r 



so that the averaged shelling function for the Fibonacci chain (and thus also for its entire LI 
class) is a(0) = 1 and a(r 2 ) = 2/ F (y*) for any non-zero distance r that is the absolute value 
of a number y E yip — A$ C Z[r]. Also, all shelling numbers a(r 2 ) are elements of Z[r], as 
can easily be seen from formula (|15l) . This has a topological interpretation, as we will briefly 
explain below for a more significant example. 

In internal space, the function / F has a piecewise linear continuation, but the function a(r 2 ) 
looks rather erratic, compare for a similar example. This is a consequence of the properties 
of the *-map, being algebraic conjugation in this case. As a mapping, it is totally discontinuous 
on L* (and also on its rational span) when the latter is given the induced topology of the 
ambient space H. In a different topology, however, this map becomes uniformly continuous, 
and it is this alternative setting, compare ^2] , which explains the appearance of the internal 
space from intrinsic data of a model set A. 

As another example, let us once more look at the circular shelling in the plane, i.e., at a 
2D model set with an open disk (radius R, centre 0) as window in 2D internal space. So, 
Q = B R (0) and, consequently, f2 — Q = B 2R (0). We have u{y) = / 2 (y*), where, due to 
rotational symmetry of the window, the function f 2 only depends on s = \y*\. Explicitly, it 
is given by 



^(*«(°> "*«('>) _ 1 arccos - 



<16) h(s) = vol(B fi (0)) = ^ arCC ° S \2R) rfV \2RJ 




Figure 1. Radial component of the frequency functions f m (s) of Eqs. (|16|) . 
(fTU|) and (|2*U)) for dimensions m = 1,2, 3, 5, 10, 50 of internal space. 
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for < s < 2R and / 2 (s) = otherwise. Fig.^contains a graph of / 2 ( s )- This function, often 
called the covariogram of the disk, is a radially symmetric positive definite function known 
as Euclid's hat, see [22l p. 100]. 

To calculate a(r 2 ), one has to sum finitely many terms of this kind, according to Eq. (|llj) . 
This situation of a 2D internal space shows up for planar model sets with n £ {5,8,12}, 
because these are the cases with <p(n) = 4. Here, one simply obtains 

(17) a(r 2 ) = c(r 2 )f 2 (s) 

where c(r 2 ) is the central shelling function of Eq. (JSJ) and s = \y*\ for any y on the shell of 
radius r. This is so because the window is a disk and the *-map sends all cyclotomic integers 
on a circle to a single circle in internal space. Consequently, the central shelling provides an 
upper bound for the average shelling in this case. 

Part of this result can be extended to arbitrary dimension. For two intersecting Tri- 
dimensional balls of radius R, the overlap consists of two congruent ball segments. The 
corresponding volume can be calculated by integrating slices (which are balls of dimension 
m — 1). Dividing by the volume of the m-ball, the covariogram becomes 

arccos( 

2V(— 4- I) r 

(18) f m (s) = r ^ ( Z + i\ / sm">)da. 

The integral can be expanded in terms of Chebyshev polynomials. For even m = 21, this 
yields 



(19) 




where Uk(x) = sin ((k + 1) arccos(x)) / sin (arccos(x)) are the Chebyshev polynomials of the 
second kind 1, Ch. 22]. For odd dimension, m = 21 + 1, one obtains the following expression 



(20) t 2e+1 (s) - 1 - 2M _ lv ^ r{e + i)jL 2kTl \£-k) 2fc+1 \2RJ 

in terms of the Chebyshev polynomials T^(x) = cos (k arccos(x)) of the first kind Ch. 22]. 
Eqs. (|18 )l — ([20 )1 are valid for < s < 2R; for distances larger than the diameter, the overlap 
vanishes, hence f m (s) = for s > 2R. Eq. (|16|) is recovered from ()19|) for £ = 1. The functions 
f m (s) for various dimensions m are shown in Fig. ^ Unfortunately, for m > 2, there is no 
simple generalisation of Eq. I|17|). because the *-map is then more complicated. 

Let us finally consider an eightfold symmetric model set in the plane, based on the clas- 
sical Ammann-Beenker or octagonal tiling, compare [21 El E] and references therein. It is 
usually described by projection from four dimensions, where we use the lattice r = \[27L A '. 
The projections ir and ir H of (|SJ) are essentially determined by compatibility with eightfold 
symmetry. In a convenient coordinatisation j^j, the images a,j, j G {1, 2, 3, 4}, of the standard 
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Figure 2. A patch of the Ammann-Beenker tiling with vertex set tIab (left) 
and the *-image of TIab inside the octagonal window in internal space (right), 
with relative scale as described in the text. 



basis vectors of the lattice have unit length in physical space, and the same is true of the cor- 
responding projections a* in internal space. Observing that Zoi + 7La<i + TLa^ + TLa^ = Z[£ g ], 
we can continue with a formulation based on the cyclotomic integers, compare |38| . For the 
Ammann-Beenker tiling, the window is then a regular octagon O of unit edge length, see 
Fig. [21 Note that the window is invariant under the symmetry group D$ of order 16. 
Explicitly, the corresponding point set in the plane is given by 



(21) 



Am 



where * is the Galois automorphism defined by £ 8 i— > £|. If we choose ^ 8 = £ = e 27 ™' 8 and 
identify K 2 with C, this gives cij = £ J 1 , 1 < j < 4, while the *-images satisfy a* = 1 ', 
compare Fig. 03 

A somewhat tedious, but elementary calculation on the basis of Fig. Ogives 
Fact 2. The covariogram of the regular octagon of edge length one is 





(22) 



/ab(*.«) 



(X-2)(x+y)x + \(l-x) + (1-y) 



(X-2)(x 2 -y 2 ) _ (A-]i,- 

4 2 
(X-2)(x-y-l)y _ x , -, 
2 2 



+ 



2 

A+2 



if X < x 
if X — y < x < X 
if l + y<x<X — y 
ifO<x<l + y 
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Figure 3. Vectors a,- in physical and a* in internal space related by the *-map. 



where A = 1 + y/2 } 



s cos a £ 



s sin a E 



0,§V2- >/2 



and where a' is the unique angle in the interval 
of the octagon. 



0, that is related to a by the D% symmetry 

□ 



A contour map of fj^{s, a) is shown in Fig.|SJ It demonstrates that the previous consider- 
ation of a circular window is actually a reasonable approximation to this case. It is sufficient 
for most applications concerning (powder) diffraction, compare |251 Ch. 3]. 

We can now calculate the averaged shelling coefficient a(r 2 ) of (jllj) explicitly for any 
distance r in /Iab- The results for all distances with < r 2 < 5 are summarised in Tabled 
They confirm the results of [Hj which had been obtained numerically. 

As an explicit example, let us consider the shortest distance in the model set. This is 
r = v2 — \/2 = 2sin(^) which is realised by the short diagonal of the rhomb. In this case, 
there are eight numbers z £ tIab — ^AB that contribute to Eq. They form a single 




Figure 4. Two overlapping regular octagons at distance s. 
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2 



1 







-1 



-2 




-2-10 1 2 



Figure 5. Contour map of the function / AB (s,a) of (j22() in internal space. 
The contours show the decrease of / AB from its maximum value / AB (0, 0) = 1 
to / AB (s,a) = outside the outer octagonal contour. 

Dg orbit; one representative is listed in Table ^ Due to the symmetry of the window, the 
contribution of each member of the orbit is the same, so it suffices to consider a representative 
and to multiply the result by the corresponding orbit length. Choosing z = 1 — £, we find 
z* = 1 — £ 3 , hence \z*\ = y2 + \/2 = 2cos(|); the corresponding angle is a = — |, hence 
a' = The distance in internal space is rather large and the overlap area correspondingly 
small. Fact |U yields x = l + y=(2 + \[2 ) /2; hence the coefficient is given by 

(23) a(2- v / 2) = 8 / AB (2 cos (f ), f ) = 8(1 - f) =4-2^2, 

compare Tabled 

The other entries in Table ^ are calculated along the same lines. Note that s can be 
calculated from r directly via s 2 = (r 2 ) , where * coincides with algebraic conjugation in 
Q(\/2 ), defined by \[2 i— > —y2. Continuing the calculations, one faces increasing complication 
with growing distance r, and, in general, one has to expect contributions from several Dg 
orbits. For r = \/3, there is still only a single orbit, this time of length 16. Hence it again 
suffices to consider a single representative, for instance z = 1 — £ — £ 2 whose *-image is 
z* = \ ((2 + \/2 ) + i (2 — y/2 )) . More generally, the standard orbit analysis reduces the sum 
in (|11|) to a formula with one contribution per Z^s-orbit, weighted with the corresponding orbit 
length. The latter is 8 whenever the corresponding angle a of Fact |2] is an integer multiple of 
^ (which corresponds to symmetry directions of the octagon), and 16 otherwise. 

The averaged shelling numbers a(r 2 ) of Table ^ as we U as the numerically determined 
values in [HJ Tab. 2], are always elements of 2Z[\/2], i.e., numbers of the form 2k + 2£\/2 
with k,£ £ Z. The formula (|13|) for the covariogram of the regular octagon only implies air 2 ) 
to be in Q(V2). However, Eqs. (jTTj) and (|T2*|) show that the averaged shelling number a(r 2 ) 
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Table 1. The averaged shelling numbers for distances r with < r 2 < 5 in 
/Iab- Representatives z are given in terms of £ = e 27n / 8 . The examples listed 
comprise a single Dg orbit each. 

r z orbit length s a' a(r 2 ) 



V2-V2 




8 


^2 + ^2 


77 

8 


4-2^2 


1 


1 


8 


1 





4 


V2 


1 + e 2 


8 


V2 





6^2-6 


Vs 




16 




arctan ( 2 ;^) 


20-12^/2 


V2 + V2 


i+£ 


8 


V2-V2 


77 

8 


36 - 22^2 


2 


2 


8 


2 





2^2-2 


V5 


2+e 2 


16 




arctan (|) 


40V2 - 56 



is expressed as a finite sum of cluster frequencies. The latter are elements of the so-called 
frequency module Fab of the Ammann-Beenker model set /Lab, i-e., the integer span of the 
frequencies of finite clusters of /Iab (i-e., of all intersections /Iab H C with C C M 2 compact). 
Since /Iab is a Delone set of finite local complexity |27[ I38j , there are only countably many 
different clusters up to translations, and only finitely many up to a given diameter. 

The frequency module has originally been calculated by means of C*-algebraic i'T-theory 
|131 119j , but can also be obtained from simpler cohomological considerations 20 . For our 
case, the result can be simplified, by a short calculation, as follows. 

Fact 3. [13| I2f)j The frequency module of the Ammann-Beenker vertex set tIab is 

F AB = I k + ^ \k,£eZ, k + £ even}. 

Thus, 8-Fab is an index 2 submodule of Z[v2]- D 

Since the window of ^4ab is a regular octagon, each cluster that contributes to the sum 
in Eq. (jllj) occurs in either 8 or 16 orientations with the same frequency. Consequently, 
the averaged shelling numbers a(r 2 ) are elements of 8-Fab, though we presently do not know 
whether they generate the full frequency module or a subset thereof. This is consistent with 
the findings of Tableland \8, Tab. 2], and establishes an interesting link between geometric 
and topological properties of model sets [HI El ■ 

A similar treatment is possible for other examples, such as the tenfold symmetric rhombic 
Penrose |8j and Tubingen triangle 6 tilings, or the twelvefold symmetric shield tiling [7]. 
However, it is clear that, for the standard model sets, the calculation becomes rather involved, 
and we presently do not know how to determine the averaged shelling function via a generating 
function such as (|TU|) in closed form, or even whether that is the most promising way to 
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proceed. For physical applications, however, often the first few terms are sufficient, and they 
can be calculated exactly from the projection method, see [HJIHIIZI for some tables and further 
results. 
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